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The  Autonomous  Ocean  Sampling  Network  Phase  Two  (AOSN-ll)  experiment  was  com  acted  in  and 
offshore  from  the  Monterey  Bay  on  the  central  California  coast  during  July  23-Septembe  6,  2003.  The 
objective  of  the  experiment  was  to  learn  how  to  apply  new  tools,  technologies,  and  analy  is  techniques 
to  adaptively  sample  the  coastal  ocean  in  a  manner  demonstrably  superior  to  traditional  it  nhodologies, 
and  to  use  the  information  gathered  to  improve  predictive  skill  for  quantities  of  interest  to  end-users. 
The  scientific  goal  was  to  study  the  upwelling/relaxation  cycle  near  an  open  coastal  bay  n  an  eastern 
boundary  current  region,  particularly  as  it  developed  and  spread  from  a  coastal  headlam  The  suite  of 
observational  tools  used  included  a  low-flying  aircraft,  a  fleet  of  underwater  gliders,  incliding  several 
under  adaptive  autonomous  control,  and  propeller-driven  AUVs  in  addition  to  moorin,  s,  ships,  and 
other  more  traditional  hardware.  The  data  were  delivered  in  real  time  and  assimilated  int  ■  the  Harvard 
Ocean  Prediction  System  (HOPS),  the  Navy  Coastal  Ocean  Model  (NCOM),  and  the  J  t  Propulsion 
Laboratory  implementation  of  the  Regional  Ocean  Modeling  System  GPL/ROMS). 

Two  upwelling  events  and  one  relaxation  event  were  sampled  during  the  experiment. '  he  upwelling 
in  both  cases  began  when  a  pool  of  cold  water  less  than  13  C  appeared  near  Cape  Ar  )  Nuevo  and 
subsequently  spread  offshore  and  southward  across  the  bay  as  the  equatorward  wind  strt  ;s  continued. 
The  primary  difference  between  the  events  was  that  the  first  event  spread  offshore  an;l  southward, 
while  the  second  event  spread  only  southward  and  not  offshore  The  difference  is  attr  buted  to  the 
position  and  strength  of  meanders  and  eddies  of  the  California  Current  System  offshore,  v  hich  blocked 
or  steered  the  cold  upwelled  water.  The  space  and  time  scales  of  the  mesoscale  variabilii  v  were  much 
shorter  than  have  been  previously  observed  in  deep-water  eddies  offshore.  Additional  pi  >cess  studies 
are  needed  to  elucidate  the  dynamics  of  the  flow 
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1.  Introduction 

The  drive  towards  integrated  coastal-ocean  observing  systems 
has  fostered  the  need  for  basic  research  not  only  on  the  individual 
elements  composing  such  systems  but  also  on  how  to  make  these 
elements  function  together  smoothly  as  a  whole.  As  a  step 
towards  providing  improved  coastal  forecasting  over  the  coming 
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decades,  the  second  Autonomous  Ocean  Sampling  Network 
(AOSN-ll)  experiment  was  conducted  in  the  Monterev  Bay  during 
late  July-early  September  2003.  The  objective  of  th  ?  Predictive 
Skill  Experiment  was  to  learn  how  to  apply  new  toi  Is,  technol 
ogies,  and  analysis  techniques  to  adaptively  sample  the  coastal 
ocean  in  a  manner  demonstrably  superior  to  tradit  onal  meth 
odologies,  and  to  use  the  information  gathered  o  improve 
predictive  skill  for  quantities  of  interest  to  end-i  sers.  These 
quantities  might  include  the  location  of  fronts  and  boundaries 
relevant  to  outfall  dispersal,  particle  trajectories  for  people 
or  pollutants  in  the  water,  and  the  biological  respor  >e  to  ocean 
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physics  in  the  form  of  bioluminescence  or  harmful  algal  blooms, 
The  goal  necessitates  an  underwater,  through-the-air,  and  over¬ 
land  networking  infrastructure  from  which  the  program  takes  its 
name. 

Although  AOSN-11  was  a  novel  collaborative  effort  involving 
state-of-the-art  autonomous  observing  platforms  and  two  data- 
driven  models  running  in  real  time  (Chao  et  al„  2009;  Haley  et  al., 
2009)  with  data  assimilation  and  ensemble  uncertainty  predic¬ 
tions  (Lermusiaux,  2007),  the  AOSN-11  project  was  influenced 
by  previous  real-time  applications  of  observing  networks  and 
modeling  systems  in  coastal  regions.  Some  of  the  early  ocean 
observing  and  prediction  system  (OOPS)  efforts  are  reviewed 
by  Robinson  et  al.  (1998)  and  Robinson  and  Lermusiaux 
(2002).  Significant  observing  and  prediction  coastal  efforts 
also  are  described  in  manuscripts  (Lynch  and  Davies,  1995; 
Mooers,  1999;  Pinardi  and  Woods,  2002)  and  review  articles 
(e.g.,  Robinson  and  Glenn,  1999;  Dickey,  2003;  Lermusiaux 
et  al.,  2006).  Without  being  exhaustive,  research  on  coastal 
OOPS  include  scientific  research  exercises  in:  Boston  Harbor  and 
Massachusetts  Bay  (Littoral  Ocean  Observing  and  Prediction 
Systems:  Robinson  and  the  LOOPS  group,  1999;  Lermusiaux, 
2001);  in  Georges  Bank  (Lynch  et  al.,  2001);  and,  in  the  Martha’s 
Vineyard  region  within  the  Coupled  Marine  Boundary  Layers  and 
Air-Sea  Interaction  Initiative  (CBLAST:  Edson  et  al.,  1999,  2006). 
AUVs  and  models  also  have  been  jointly  utilized  within  sustained 
efforts  in  the  Gulf  of  Maine  such  as  the  Gulf  of  Maine  Ocean 
Observing  System  (GoMOOS,  Bogden  et  al.,  2001 )  and  the  Regional 
Association  for  Research  on  the  Gulf  of  Maine  (RARGOM:  Runge, 
2005).  Several  of  these  research  activities  in  the  Northeast 
involved  coupled  physical-biological  observing  and  modeling 
research,  from  determining  dominant  scales  to  targeted  process 
studies,  e.g..  Ecology  of  Harmful  Algal  Blooms/Global  Ocean 
Ecosystem  Dynamics  (ECOHAB/GLOBEC:  Anon.,  2002.  2005; 
Wiebe  et  al.,  2001)  and  LOOPS  (Besiktepe  et  al.,  2003),  but  for 
biological-physical  data  collection  and  modeling  in  real  time, 
interesting  challenges  remain  (e.g.,  deYoung  et  al ,  2004 ).  The  first 
Long-term  Ecosystem  Observatory  (LEO)  research  with  advanced 
ocean  and  observing  prediction  systems  has  been  maintained  off 
the  coast  of  New  Jersey,  with  numerous  data-model  integrated 
coastal  studies  and  results  (e.g.,  Schofield  et  al.,  2002;  Wilkin 
et  al.,  2005).  Important  research  has  also  being  carried  out  in  the 
Gulf  or  Mexico  and  its  larger  region  of  influence  (Mooers  et  al., 
2005).  Numerous  scientific  results  also  have  been  obtained  off  the 
coast  of  Oregon  as  a  result  of  coupled  observation  and  modeling 
experiments  (Barth  and  Wheeler,  2005).  All  of  these  research 
efforts  are  now  better  coordinated  and  developed  at  the  national 
level,  within  programs  such  as  the  Integrated  and  Sustained  Ocean 
Observing  Systems  (Ocean.US,  2002,  2004).  At  the  international 
level,  several  exercises  were  sponsored  in  multiple  regions  by  the 
North  Atlantic  Treaty  Organization  within  the  Rapid  Environ¬ 
mental  Assessments  and  Rapid  Response  exercises  (Bovio  et  al., 
2003;  Robinson  et  al.,  2003;  Onken  et  al.,  2002,  2005).  Specifically, 
the  first  adaptive  sampling  based  on  real  time  data-driven 
uncertainty  predictions  was  carried  out  within  in  a  NATO  research 
exercise  in  the  Strait  of  Sicily  region  (Robinson  et  al.,  1999, 
Lermusiaux,  1999). 

While  many  of  the  objectives  of  the  AOSN-11  experiment  were 
practical  in  nature,  the  AOSN-11  science  team  formulated  a  list  of 
scientific  goals  as  well  including: 

•  Observe  and  predict  the  development  and  movement  of 
upwelling  fronts  and  ocean  eddies. 

•  Determine  the  importance  of  the  wind-stress  curl  in  driving 
upwelling  and  the  three-dimensional  circulation  in  the 
Monterey  Bay. 


•  Study  the  relationship  between  micro-scale  atmospheric  jets 
and  the  formation  of  oceanic  cold  plumes  off  coastal  headlands 
and  promontories. 

•  Examine  the  nutrient  supply  processes  in  the  Monterey  Bay 
including  upwelling,  advection,  and  mixing. 

•  Evaluate  and  improve  numerical  models  of  both  the  atmo 
sphere  and  the  coastal  ocean. 


The  Monterey  Bay  provides  an  ideal  location  for  such  an 
experiment  due  to  its  close  proximity  to  several  large  research 
institutions,  ship  availability,  and  the  wealth  of  environments  and 
scientific  problems  presented.  The  wide  continental  shelves 
within  and  to  the  north  of  the  bay  provide  a  strong  contrast  to 
the  nearby  deep  water  within  the  Monterey  Lay  Submarine 
Canyon  (MSC)  and  the  narrower,  steeper  shelf  to  the  south. 
The  two  upwelling  centers  at  Point  Sur  and  Cape  Ano  Nuevo  are 
both  within  a  day’s  steam  by  research  vessel.  Much  of  the 
communications  infrastructure  including  underse  i  cables,  semi 
permanent  moored  surface  buoys,  and  a  wireles  Internet  net 
work  was  already  established  by  the  Monterey  Bay  Aquarium 
Research  Institute  and  various  local  Navy  entitie5  A  network  of 
coastal  HF  radars  was  already  in  place  momtori  ig  the  surface 
currents  in  real  time.  The  installation  of  the  addn  onal  elements 
required  for  the  AOSN-11  experiment  to  the  exiving  observing 
system  in  the  bay  was  straightforward. 

The  mesoscale  variability  near  the  Montere  r  Bay  can  be 
succinctly  described  as  the  interplay  between  the  upwelling 
centers  rooted  at  headlands  to  the  north  and  so  ith  of  the  bay 
and  a  persistent,  anti-cyclonic  California  Current  meander  just 
offshore  of  the  bay  itself  (Ramp  et  al.,  2005).  The  upwelling 
centers  are  clearly  driven  by  the  prevailing  north\  esterly  winds, 
which  are  strongest  during  March-June  (Nelson,  177).  The  flow 
off  Ano  Nuevo  is  sometimes  bifurcated  (Rosenfel  1  et  al.,  1994), 
with  some  of  the  cold,  salty  upwelled  water  flowii  g  offshore  and 
some  flowing  south  over  the  Monterey  Canyor  This  pattern 
seems  to  depend  on  the  wind-stress  duration:  1  he  cold  water 
spreads  southward  only  after  a  strong,  steady  nortl  westerly  wind 
has  been  blowing  for  about  6-7  days  (Ramp  et  al.  2005). 

To  the  south  of  the  bay,  a  strong  upwelling  ce  iter  located  at 
Point  Sur  generates  cold  filaments  moving  offshoi  1  in  the  upper 
100m  (Breaker  and  Mooers,  1986;  Tisch  et  al.,  1  92;  Rosenfeld 
et  al.,  1994;  Ramp  et  al.,  1997).  On  the  coi  tinental  shelf 
there,  currents  were  equatorward  during  north\  esterly  winds 
and  poleward  during  wind  reversals  and  relaxatii  ns  (Ramp  and 
Abbott,  1998).  These  authors  found  that  the  veiical  structure 
of  equatorward  currents  was  consistent  with  the  s  perposition  of 
the  local  wind  forced  Ekman  spiral  and  the  ale  igshore  geos- 
trophic  flow  due  to  the  set-down  at  the  coast.  The  oleward  flows 
during  relaxations  were  dynamically  consistent  w  th  a  poleward 
alongshore  pressure  gradient  force,  as  observed  at  >ther  locations 
off  California  (Winant  et  al.,  1987;  Lentz,  1987).  The  ,e  currents  are 
consistent  with  the  view  that  cold  surface  water  i  i  the  center  of 
the  Monterey  Bay  is  advected  from  the  north  of  tl  e  bay  and  not 
the  south. 

The  anti-cyclonic  California  Current  (CC)  meant  er,  also  some¬ 
times  referred  to  as  the  Monterey  Bay  Eddy  (MBE),  is  a  frequently 
observed  feature  of  the  region  during  the  up\  'elling  season 
(Rosenfeld  et  al.,  1994;  Paduan  and  Rosenfeld,  199  i;  Ramp  et  al., 
2005).  Drifters  suggest  that  this  feature  is  assoc  ated  with  the 
larger-scale,  meandering  California  Current  System  and  is  not 
locally  generated  (Brink  et  al.,  1991).  Moored  obs<  rvations  show 
that  the  feature  is  deep,  with  coherent  flows  greatt  r  than  20  cm/s 
exceeding  1000  m  depth  (Ramp  et  al.,  1997).  The  feature  is  warm 
and  fresh  relative  to  the  local  waters,  displaying  t  ie  influence  of 
the  Pacific  Subarctic  Water  (Rosenfeld  ct  al.,  1994).  The  MBE  is 
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thus  an  oceanic  feature  and  is  clearly  not  generated  by  the  local 
wind  stress,  although  it  does  respond  to  it.  Time  series  of  sea- 
surface  temperature  from  aircraft  and  AVHRR  images  show  that 
the  MBE  moves  rapidly  onshore  during  wind  relaxation  events 
and  retreats  back  offshore  when  the  winds  reintensify  (Rosenfeld 
et  al.  1994).  Moored  observations,  however,  also  have  shown  an 
onshore  and  southward  translation  of  the  MBE  without  a  wind 
relaxation  (Ramp  et  al.,  1997).  The  mechanism  for  this  across- 
shore  eddy  translation  is  thus  not  fully  understood. 

Coursing  poleward  beneath  the  upwelling  centers  and  CC 
meander  lies  the  sub-surface  California  Undercurrent  (CUC).  The 
CUC  is  observed  over  the  continental  slope  all  along  the  west  coast 
of  the  United  States  (Pierce  et  al.,  2000)  and  is  a  ubiquitous  feature 
of  all  moored  data  sets  both  north  and  south  of  the  Monterey  Bay. 
To  the  north  off  the  Farallon  Islands,  the  undercurrent  was  present 
most  of  the  year  but  was  not  coherent  with  local  wind  forcing 
(Noble  and  Ramp,  2000).  To  the  south  off  Point  Sur,  the  currents 
below  100  m  depth  are  more  often  poleward  than  equatorward 
and  frequently  exceed  30 cm/s  at  100  m  (Chelton,  1984;  Wickham 
et  al.,  1987;  Tisch  et  al.,  1992;  Ramp  et  al..  1997).  The  poleward 
flow  off  Point  Sur  is  pulse-like  at  very  low  frequencies  (3-4 
months)  (Ramp  et  al.,  1997).  How  this  pulse-like  flow  relates 
to  cross-shore  translations  of  the  MBE  and  potential  “blocking” 
of  the  CUC  by  the  eddy  is  of  interest  but  poorly  understood. 

The  circulation  in  the  bay  itself  is  usually  cyclonic  (Paduan  and 
Rosenfeld.  1996;  Paduan  and  Cook.  1997;  Ramp  et  al.,  2005) 
during  both  upwelling  and  downwelling  conditions.  During 
upwelling  events,  very  warm  (>16  C)  water  is  observed  in  the 
northeast  corner  of  the  bay,  in  the  wind  shadow  behind  the  Santa 
Cruz  mountains  (Graham.  1993;  Graham  and  Largier,  1997;  Ramp 


et  al.,  2005).  Model  salinity  fields  at  the  surface  and  100-m  levels 
show  that  the  bay  fills  with  higher-salinity  water  stemming  from 
the  Point  Ano  Nuevo  upwelling  center  under  these  conditions. 
During  downwelling,  the  cyclonic  circulation  was  more  barotropic 
over  the  continental  shelf  than  during  upwelling  events.  The 
source  water  for  both  the  surface  and  100-m  levels  was  the  colder, 
fresher  California  Current  water  offshore,  which  had  advected 
southward  well  past  Point  Pinos  during  the  previous  upwelling 
event  and  entered  the  bay  from  the  south  (Ramp  et  al..  2005) 


2.  Data  and  methods 

An  impressive  array  of  sampling  platforms  was  assembled 
for  AOSN-II  including  ships,  aircraft,  gliders,  and  pro  )eller-driven 
AUVs,  with  each  deployed  to  optimize  its  sampli  ig  strengths 
(Fig.  1 ).  Gliders  are  buoyancy-driven  vehicles  that  are  slow  but  can 
remain  deployed  for  several  months  at  a  time  (Belli  gham  et  al., 
2009).  Propeller-driven  vehicles  can  sustain  speeds  o  6  knots  but 
only  for  14  h  or  less  per  deployment.  Ships  collect  he  highest- 
accuracy  data  and  can  collect  water  from  througho  it  the  water 
column  for  chemical  analyses.  Aircraft  can  obtain  a  “  .napshot”  of 
the  sea-surface  conditions  in  5  h  or  less,  but  with  1  ss  accuracy 
than  ships. 

Five  Spray  gliders  (Sherman  et  al.,  2001;  Davis  t  al.,  2002) 
were  deployed  during  AOSN  II  from  July  23  to  Aug  st  24.  2003 
along  5  linear  tracks  covering  from  Point  Ano  Nuevo  in  the  north 
to  Point  Sur  in  the  south,  out  to  roughly  100km  offs  lore  (Fig.  1 ). 
Beyond  instrument  development  and  testing  operatn  nal  capabil 
ities,  the  goal  for  these  gliders  was  to  contribi  te  data  for 
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Fig.  1.  Sampling  locations  for  autonomous  vehicles  (left)  and  other  assets  such  as  ships,  aircraft,  and  HF  radars  (right)  during  AOSN-II.  Green  tracks  were  ob  lined  from  the 
Spray  vehicles  and  the  white  tracks  from  the  Slocums.  The  blue  ‘V*  track  was  occupied  at  night  by  the  CalPoiy  REMUS  vehicle.  The  other  red  tracks  were  (  cupied  by  the 
MBARI  DORADO  vehicle  and  the  NPS  REMUS  vehicle  The  dark  blue  track  on  the  right  is  the  aircraft  flight  path,  and  the  thin  dotted  line  with  station  locatio  s  is  the  POINT 
SUR  ship  survey.  The  vectors  in  the  Monterey  Bay  are  the  surface  currents  as  sensed  by  coastal  HF  radars.  Both  sampling  schemes  are  overlaid  on  satellite  AVI  RR  sea  surface 
temperature  (SST),  as  indicated  by  the  color  bar. 


S.R,  Ramp  et  al.  (  Deep-Sea  Research  II  56  (2009)  68  86 


71 


assimilation  and  synoptic  mapping  and  to  define  the  scales  and 
patterns  of  variability.  These  gliders  collected  a  total  of  2075 
vertical  profiles  to  400  m  depth  (or  the  bottom)  with  a  nominal 
spatial  resolution  of  2.3  km.  The  resolution  is  a  function  of  bottom 
depth  and  higher  resolution  was  obtained  in  shallower  water.  In 
addition  to  temperature,  salinity,  and  water  clarity,  each  glider 
also  returned  an  estimate  of  the  vertically  averaged  horizontal 
current  during  each  profile,  obtained  as  the  difference  between 
actual  GPS  position  and  dead  reckoning.  A  Spray  glider  can 
typically  make  only  about  22-25 cm/s  of  forward  progress 
through  the  water,  which  necessitated  the  across-shore  racetracks 
to  avoid  fighting  the  primarily  alongshore  current  in  the  region. 
It  took  each  glider  nominally  4  days  to  transit  one  way  along  each 
100-km  racetrack. 

Overall  Spray  glider  performance  was  excellent  during  AOSN  II. 
The  gliders,  once  deployed  required  no  servicing  and  lasted 
for  the  duration  of  the  experiment  on  one  set  of  batteries.  The 
limiting  factor  on  good  data  collection  turned  out  to  be  biofouling 
rather  than  battery  life.  For  additional  details,  see  the  accom¬ 
panying  technical  article  (Bellingham  et  al.,  2009). 

Ten  Slocum  gliders  profiled  to  200m  depth,  and  were  thus 
deployed  closer  to  shore  than  the  Sprays  (Fig.  1).  Four  Slocum 
gliders  ran  closed  boxes  (or  ’cells”)  rather  than  linear  tracks,  and 
were  equipped  with  fluorescence  and  photosynthetically  available 
radiation  (PAR)  sensors  in  addition  to  the  variables  sensed  by  the 
Spray.  A  total  of  11,341  profiles  yielding  nominally  800  m  along  - 
track  resolution  were  obtained  over  the  life  of  the  experiment. 
The  endurance  of  the  Slocums  was  11-14  days,  about  half 
that  expected,  and  some  recovery  and  re-deployment  was  thus 
required  to  cover  the  lifetime  of  the  experiment.  The  Slocums  also 
experience  problems  with  biofouling,  but  the  shorter  deploy¬ 
ments  allowed  for  cleaning  of  the  sensors  to  minimize  this 
problem 

One  of  the  most  innovative  aspects  of  the  AOSN-II  experiment 
was  adaptive  sampling  of  ocean  features  using  gliders  (Leonard 
et  al.,  2007).  Three  to  four  Slocum  gliders  were  used  for  this 
purpose.  Three  demonstrations  of  autonomous  adaptive  control 
were  successfully  completed,  one  in  the  northern  part  of  the  bay 
across  an  upwelling  front,  the  second  in  the  southern  part  of  the 
bay,  and  the  third  following  a  surface  drifter  (Fiorelli  et  al.,  2006). 
For  the  first  experiment  during  August  6-7,  2003,  three  gliders 
maintained  a  triangular  pattern  3  km  on  a  side  while  propagating 
NW  along  a  linear  path.  The  triangular  formation  was  chosen  so 
that  m-situ  estimates  of  gradients  (e.g.,  in  temperature)  along  the 
path  of  the  center  of  the  formation  could  be  computed  in  near-real 
time  from  the  gliders’  scalar  (temperature)  measurements.  The 
gliders  could  be  programmed  to  follow  the  gradient  if  desired. 
The  3-km  sensing  array  resolution  was  designed  with  respect 
to  previously  computed  scales  in  Monterey  Bay.  The  formation 
error  was  423  m,  or  14%  of  the  desired  spacing  with  a  standard 
deviation  of  159  m  For  the  second  demonstration  on  August 
16-17,  2003,  some  additional  difficulty  was  added.  The  triangle 
formation  was  asked  to  travel  a  zigzag  path  across  an  upwelling 
front,  and  to  shrink  the  size  of  the  formation  from  6  to  3  km  on  a 
side  enroute.  The  objective  was  to  demonstrate  the  ability  of  the 
glider  control  system  to  adjust  the  formation  size,  i.e.,  the  sensing 
array  resolution,  in  response  to  changing  scales  in  the  sampled 
field(s)  and  to  provide  gradient  estimates  at  different  resolutions. 
This  exercise  was  also  a  success:  The  glider  formation  error  was  a 
function  of  the  size  of  the  triangle,  about  394  m  for  the  6-km 
triangle  and  270  m  for  the  3  km  one.  Finally,  on  August  23.  2003,  a 
single  glider  was  asked  to  follow  a  surface  drifter  in  real  time, 
making  zigzags  across  its  projected  path.  The  goal  was  to 
demonstrate  that  a  glider  or  glider  formation  could  collect  scalar 
measurements  (or  gradient  estimates)  both  across  and  below 
tracer  paths.  This  demo  was  qualitatively  a  success,  but  there  were 


some  problems  based  on  the  different  velocity  fields  sensed  by  the 

two  instruments:  The  drifter  measures  the  surface  velocity  while 
the  glider  measures  the  current  vertically  integrated  over  the 
up/down  path,  making  coincident  path  prediction  difficult.  There 
are  other  (easier)  ways  to  do  this,  such  as  having  the  glider  follow 
the  drifter  with  a  two-hour  lag.  which  may  be  more  useful  in 
different  applications. 

Propeller-driven  AUVs  such  as  the  commerrially  available 
REMUS,  the  MBARI  manufactured  DORADO,  and  the  NPS  in  house 
vehicle  ARIES  typically  have  about  14  h  duration  at  4-6  knots. 
The  vehicles  also  have  a  much  larger  payload  than  gliders  and  can 
thus  carry  more  sophisticated  sensors.  The  larg  *st  number  of 
transects  were  run  by  the  California  Polytechnic  5 rate  University 
(CalPoly)  REMUS  vehicle  along  two  transects  in  th  northern  half 
of  the  Monterey  Bay  The  goal  of  sensing  biolumine  cence  dictated 
that  these  REMUS  runs  be  done  in  the  dark,  and  the  efore  9  seven 
hour  missions  were  completed  along  45  km  trar  sects  between 
9  pm  and  4  am  each  night  (blue  “V"  in  Fig.  1).  Th  ese  operations 
began  near  Santa  Cruz,  ran  out  to  MBARI  surface  buoy  Ml,  then 
returned  to  the  NW  corner  of  the  bay  (Molini  et  al.,  2005; 
Shulman  et  al.,  2005).  The  vehicle  followed  a  saw  t  >oth  pattern  to 
40m  depth  sampling  with  a  CTD,  transmissometi  r,  fluorometer, 
and  ADCP  in  addition  to  the  bioluminescence.  Four  transects  were 
also  run  from  Monterey  to  Ml  using  the  NPS  F  EMUS  vehicle, 
sampling  temperature  only.  The  NPS  ARIES  veh  cle  was  used 
exclusively  to  fetch  data  from  a  bottom-mountt  i  ADCP  using 
acoustic  modems  and  did  not  collect  environment  il  data. 

The  largest  AUV  used  during  AOSN-ll  was  the  IV  BARI  DORADO 
AUV.  This  vehicle  ran  along  a  line  between  Moss  Landing 
and  MBARI  buoy  M2,  profiling  to  250  m  depth  (Fig.  ).  Instruments 
on  board  included  a  CTD,  fluorometer,  oxygen  and  itrate  sensors, 
bioluminescence,  and  ADCP.  The  run  was  madt  operationally 
every  3  weeks,  but  the  frequency  was  increased  luring  August 
2003  to  1  section  per  week.  Two  additional  DOF  \DO  transects 
also  were  conducted  offshore.  A  clear  advantage  )f  this  vehicle 
was  the  sampling  depth,  which  was  very  helpful  or  model  data 
assimilation.  Its  large  size  requires  a  specializ'  d  vehicle  for 
recovery. 

The  airborne  measurements  during  AOSN-1!  vere  obtained 
using  the  TWIN  OTTER  aircraft  owned  and  operated  by  the  Center 
for  Interdisciplinary  Remotely  Piloted  Aircraft  St  dies  (CIRPAS) 
at  the  Naval  Postgraduate  School.  Fourteen  flights  v  ere  conducted 
between  August  4  and  September  5,  2003  at  n  altitude  of 
nominally  33  m  above  the  sea  surface.  The  lov  altitude  was 
chosen  to  pass  beneath  the  stratus  cloud  deck  th,  t  often  covers 
the  Monterey  Bay  in  summer,  and  to  minimi2  *  the  vertical 
divergence  in  the  heat  and  momentum  flux  observations.  The 
flight  path  (Fig.  1,  blue  line)  covered  all  of  the  >ay  and  some 
distance  to  the  north  and  south  out  to  100  km  offshore.  Two 
vertical  saw-tooth  paths  to  1500  m  also  were  co  iducted  along 
two  slightly  longer  transects  at  the  northern  end  of  the  flight  path 
and  along  the  center  of  the  bay  (Fig.  1).  These  )aths  were  to 
sample  the  three-dimensional  structure  of  the  narine  atmo¬ 
spheric  boundary  layer  over  the  bay  (Kalogiros  et  a  .,  2006,  Wang 
et  al.,  2007).  The  atmospheric  parameters  sampled  by  the  aircraft 
included  air  temperature  and  dew  point,  atmosp  eric  pressure, 
wind  speed  and  direction,  the  turbulent  fluxe<  of  heat  and 
momentum,  and  total  particle  number,  as  well  «  s  aerosol  and 
cloud/fog  droplet  size  distributions.  In  the  oce^  i,  sea-surface 
temperature  (SST),  ocean-leaving  radiance  at  19  wavelengths, 
sun  glint,  and  surface  roughness  were  observed  The  SST,  air 
temperature  and  dew  point,  and  wind  data  were  made  available 
for  model  assimilation  in  near-real  time,  wh  le  the  other 
parameters  required  additional  post-processing. 

Coastal  HF  radars  for  mapping  ocean  surface  currents  have 
been  proliferating  along  the  central  coast  for  a  p  riod  of  years. 
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starting  with  three  standard-range  (13  MHz)  systems  in  the 
Monterey  Bay  located  at  Point  Pinos,  Moss  Landing,  and  Santa 
Cruz.  The  network  has  recently  been  extended  to  the  north  with 
long-range  (5  MHz)  systems  at  Big  Creek,  Pescadero,  and  Montara. 
All  the  bay  systems  were  operational  during  AOSN-U  and  provided 
surface  currents  in  2-3  km  bins  for  model  data  assimilation 
(Paduan  and  Shulman,  2004;  Shulman  and  Paduan,  2009)  and 
Lagrangian  Coherent  Structures  (LCS)  computations. 


heat  fluxes,  and  assimilated  the  aircraft  SST,  glider  T  and  S.  and 
AUV  data.  The  JPL/ROMS  model  additionally  assimilated  T  and  5 
data  from  MBAR1  moorings  Ml  and  M2.  None  of  the  models 
assimilated  velocity  from  moorings  or  HF  radar.  A  table  with 
additional  details  of  the  three  models,  the  data  assimilation 
schemes  used,  and  how  they  were  forced  may  be  found  on 
the  web  at:  (http://www.princeton.edu/~dcsl/asap/ASAP_ROMS 
HOPS_NCOMICON_summary_MB06.doc). 


3.  The  numerical  models 


4.  Results 


A  key  aspect  of  AOSN-II  program  was  to  assimilate  the  data 
collected  in  near-real  time  into  mesoscale  ocean  models  to 
forecast  conditions  for  the  following  day,  and  adapt  the  sampling 
scheme  accordingly.  In  addition  to  real-time  operations,  the  models 
also  were  used  to  help  elucidate  the  dynamics  of  the  flow.  Three 
models  were  used  for  these  purposes,  namely  the  Harvard  Ocean 
Prediction  System  (HOPS)  (Robinson.  1999),  the  jet  Propulsion 
Laboratory  implementation  of  the  Regional  Ocean  Modeling 
System  (JPL/ROMS)  (Shchepetkin  and  McWilliams,  2004),  and 
the  Navy  Coastal  Ocean  Model/Innovative  Coastal  Ocean  Obser¬ 
ving  Network  (NCOM/ICON)  (Shulman  et  al.,  2002,  2004).  While 
conceptually  similar  in  scale  (1. 5-2,5  km  in  the  horizontal  x  32 
sigma  levels  in  the  vertical)  and  operation,  the  models  differed 
significantly  in  terms  of  the  mixing  schemes  employed,  initial 
conditions  used,  spin-up  time,  and  open  boundary  conditions 
provided  by  the  larger-scale  models  they  were  nested  within.  All 
three  models  were  forced  using  the  Coupled  Air  Ocean  Modeling 
and  Prediction  System  (COAMPStM)  surface  wind  stresses  and 


4.1.  Atmospheric  forcing 

The  weather  presented  an  interesting  mix  of  conditions  that 
were  well  suited  to  the  experimental  objectives  di  ring  August 
2003.  Time  series  from  MBAR1  buoy  M2,  located  at  36  42  N. 
122  23  W  at  the  mouth  of  the  bay,  show  that  the  fir?  t  two  flights 
during  August  4-5  took  place  during  a  wind  rela  ation  event 
when  winds  were  weakly  poleward  at  less  than  5  m/ .  (Fig.  2).  An 
extended  period  of  strong  (commonly  10  m/s),  stead' ,  upwelling 
favorable  northwesterly  winds  followed  during  A  jgust  6-19, 
which  were  quite  typical  for  this  time  of  year.  Five  WIN  OTTER 
flights  sampled  this  first  upwelling  event.  Another  po  ?ward  event 
took  place  during  August  20-21,  followed  by  a  secoi  d  upwelling 
event  from  August  23-31.  Three  flights  were  executed  during  the 
relaxation  and  two  during  the  second  upwelling  even  September 
1-3  was  very  calm  and  warm  in  the  Monterey  E  ly  and  was 
characterized  by  nearly  flat  seas  and  strong  surface  eating.  This 
time  period  was  unfortunately  not  sampled  by  the  air  raft.  During 
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Fig.  2.  Low-pass  filtered  surface  wind  vectors  from  Buoy  M2,  located  at  36  42'N.  122  ’23  W.  maintained  and  operated  by  the  Monterey  Bay  Aquarium  Re  earch  Institute 
(MBAR1).  The  gray  vertical  bars  indicate  the  time  of  the  NPS  TWIN  OTTER  aircraft  over-flights. 
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Aircraft  Wind  COAMPS  3km  Wind 
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Aircraft  Curl  COAMPS  3km  Curl 


Fig.  3.  A  comparison  of  th wind,  wind  stress,  and  the  wind-stress  curl  between  the  observed  values  sampled  by  the  aircraft  and  the  output  from  the  3-km  (DAMPS1  M 
model  on  August  15,  2003  during  an  upweiling  event.  Common  scales  have  been  adopted  for  both  pictures.  The  wind  speed  displayed  as  vectors  are  ove  laid  on  the  color- 
coded  wind  stress  as  indicated  by  the  left  side  of  the  color  bar.  The  wind-stress  curl  translated  into  upweiling  velocities  in  m/d  are  indicated  by  the  lower  anel  dots  and  the 
right-hand  scale  of  the  color  bar  The  black  circle  in  the  lower  panels  indicates  the  size  of  the  domain  used  to  compute  the  curl  of  the  wind  stress. 


September  4-6,  weakly  upwelling-favorable  winds  returned  to 
the  area  once  again,  and  the  final  two  aircraft  over-flights  were 
conducted. 

Since  the  spatial  scales  of  the  atmospheric  forcing  are  large 
compared  to  the  size  of  the  Monterey  Bay  (22  x  40  km),  the  winds 
observed  at  buoy  M2  are  usually  representative  of  the  offshore 
winds  elsewhere  along  the  coast  (Fig.  3).  Within  the  bay  itself, 
however,  the  surrounding  mountains  create  high  spatial  varia¬ 
bility  in  the  surface  wind  stress.  Of  particular  interest  are  the 
weak  winds  behind  the  mountains  in  the  NE  corner  of  the  bay, 
and  a  small-scale  atmospheric  jet  coming  through  a  gap  in  the 
mountains  off  Santa  Cruz  (Fig.  3).  This  figure  from  August  15  was 
sampled  9  days  into  the  upweiling  event  and  was  typical  of  the 
others  during  this  time.  The  aircraft  and  the  Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  System  (COAMPS)  (Doyle  et  al.. 
2009)  both  show  the  upwelling-favorable  winds  turning  slightly 
into  the  bay.  Both  also  show  the  weak  (blue)  winds  in  the  NE 
corner  of  the  bay.  The  wind  jet  off  Santa  Cruz  was  slightly  clearer 
in  the  model  output  than  the  observational  data,  as  was  the  strong 
positive  curl  on  the  inshore  side  of  the  jet  (Fig,  3,  bottom). 
The  location  and  scale  of  this  feature  closely  resembled  a  similar 
one  observed  on  August  17,  2000  (Ramp  et  al.,  2005).  While  that 
jet  was  noticeably  hot  and  dry  (their  Fig.  4)  this  was  less  obvious 
in  the  2003  data.  This  is  perhaps  due  to  the  higher  altitude 
(133  m)  of  the  year  2000  flights.  The  2003  atmospheric  tempera¬ 
tures  at  33  m  (not  shown)  more  closely  followed  the  sea-surface 
temperature  and  were  cool  rather  than  warm.  There  was  a 


minimum  in  the  dew  point  temperature  just  o  f  Santa  Cruz, 
however,  indicating  slightly  dryer  conditions  th  re.  A  typical 
poleward  wind  event  from  August  20,  2003  (Fig.  4) « lowed  overall 
weaker  winds,  weaker  gradients,  and  less  wind-str  ss  curl  than  a 
typical  equatorward  wind  event.  The  wind  stress  w  is  still  weaker 
in  the  NE  corner  of  the  bay  than  elsewhere,  and  th  •  wind  turned 
into  the  bay  from  the  south  rather  than  the  north. 


4.2.  Development  and  relaxation  of  the  upweiling  p!  me 

A  time  series  of  12  realizations  of  SST  and  wind  is  sampled  by 
the  aircraft  shows  the  ocean’s  mesoscale  res  onse  to  the 
equatorward  and  poleward  wind-stress  events  d;  ring  summer 
2003  (Fig.  5A-L).  August  4  was  the  first  day  of  a  re  axation  event 
following  a  weak,  2-day  upweiling  event.  The  win  is  were  weak 
and  disorganized  and  some  traces  of  coastal  upwe  ing  were  still 
visible  ( <  14  C)  off  Point  Ano  Nuevo.  The  warmest  v  ater  exceeded 
17  C  offshore  and  in  the  southern  portion  of  the  Me  iterey  Bay.  By 
the  following  afternoon  (Fig.  5B),  the  entire  se  i  surface  had 
warmed  by  about  2  C  under  weak,  slightly  onsho  e  winds,  with 
the  exception  of  off  Point  Sur  where  the  SST  remained  about  the 
same.  There  was  some  indication  of  locally  upwe  ling-favorable 
winds  off  Point  Sur,  which  perhaps  balanced  the  local  surface 
heating. 

The  next  five  panels  (Fig.  5C-G)  show  the  devel  )pment  of  the 
largest  upweiling  event.  The  August  6  image  was  obtained  on 
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Aircraft  Cud  COAMPS  3km  Curl 

Fig,  4.  As  in  Fig.  3  except  on  August  20,  2003  during  a  relaxation  event 


the  first  day  of  strong  (10  m/s)  upwelling-favorable  winds  The 
upwelling  was  just  starting  (yellow  contour  <15  )  at  both  Point 
Ano  Nuevo  and  Point  Sur.  Note  that  the  aircraft  survey  was 
designed  to  focus  on  the  Ano  Nuevo  region  and  most  of  the  Point 
Sur  upwelling  plume  was  located  to  the  south  of  the  flight  path 
and  was  not  observed.  The  discussion  hereafter  focuses  on  the 
growth  and  spreading  of  the  Ano  Nuevo  plume.  Four  days  later 
(Fig.  5D)  the  upwelling  center  had  enlarged  and  cooled  consider¬ 
ably  to  less  than  13  C.  There  was  some  evidence  of  plume 
bifurcation  (Rosenfeld  et  al..  1994)  with  some  of  the  14  C  water 
heading  offshore  and  some  heading  south.  By  the  next  day 
however  (Fig.  5E),  southward  advection  seemed  to  be  winning 
with  all  the  13  and  14  C  water  heading  south  and  connecting  up 
with  the  cold  water  from  the  Point  Sur  upwelling  center.  This  lead 
to  a  detached  parcel  of  13  C  water,  which  advected  south  on 
August  13  (Fig.  5F).  Note  that  the  CC  meander,  as  indicated  by  the 
17  C  contour,  remained  well  offshore  during  this  entire  event, 
which  allowed  the  cold  water  to  follow  an  unobstructed  path 
towards  the  south.  On  August  15,  following  9  days  of  continuous 
upwelling,  the  winds  were  blowing  at  15  m/s  and  the  coldest 
water  (<12  C)  was  observed  off  both  Ano  Nuevo  and  Point  Sur. 
The  cold  water  from  the  two  centers  merged  together  and  formed 
a  continuous  band  of  minimum  SST  roughly  along  122  06  W  that 
separated  the  warm  water  offshore  from  the  also  warm  water  in 
the  Monterey  Bay. 

The  2-week  upwelling  event  was  followed  by  3  days 
of  poleward  winds  at  about  8  m/s  on  August  20-22  (Figs.  2 
and  5H-J).  The  ocean’s  adjustment  to  the  change  in  surface  forcing 
was  remarkably  rapid  and  shows  two  noteworthy  features 
(Fig.  5H).  First,  the  boundary  between  the  warm  California 
Current  water  moved  onshore  to  within  just  a  few  kilometers 
of  the  coast  near  Davenport,  just  north  of  Santa  Cruz,  CA.  This 


dramatically  reduced  the  amount  of  cold  water  at  thi  surface  in 
the  northern  portion  of  the  domain.  Second,  a  cc  d  ’’squirt” 
formed  in  the  southern  half  of  the  domain,  centered  near  Point 
Pinos  in  Monterey.  Since  the  wind  generally  reverses  rom  south 
to  north  along  the  coast,  we  speculate  that  this  feati  re  resulted 
from  convergence  between  poleward  currents  along  tl  e  southern 
Big  Sur  coast  and  currents  that  were  still  equatorward  3  the  north 
of  this  region.  This  feature  was  still  observable  as  a  It  ns  of  15  C 
water  55  km  offshore  on  August  21  (Fig.  51)  but  wai  ultimately 
pushed  back  onshore  by  the  continuing  incursion  of  (  C  water  on 
the  22nd  (Fig.  5J). 

The  final  two  panels  of  Fig.  5  show  the  return  tc  upwelling 
conditions,  as  represented  by  two  realizations  on  Auj  ust  25  and 
29  (Fig.  5K,  L).  The  25th  was  just  2  days  into  the  ev<  nt  and  not 
surprisingly  resembles  August  6,  with  weak  upwelling  ind  <15  C 
water  off  Cape  Ano  Nuevo.  August  29,  6  days  into  the  event  with 
winds  >15  m/s  closely  resembles  August  15.  One  significant 
difference  would  be  that  the  CC  was  closer  to  shore  <  n  the  29th 
than  the  15th,  which  confined  the  cold  water  closer  t  shore  and 
resulted  in  strictly  southward  rather  than  offshore  spreading. 
These  differences  in  the  two  upwelling  events  were  als  evident  in 
the  output  from  the  NCOM/1CON  model  (Fig.  6).  The  bCOM-ICON 
model  is  based  on  the  Navy  Coastal  Ocean  Model,  a  id  is  triply 
nested  inside  of  the  global  and  regional  (Californ  a  Current) 
NCOM-based  models  (Shulman  et  al.,  2004).  The  mo<  el  is  called 
NCOM  ICON  due  to  the  fact  that  initial  development  o  the  model 
started  under  the  National  Oceanic  Partnership  Progi  im  (NOPP) 
Innovative  Coastal-Ocean  Observing  Network  (ICC  N)  project 
(Ramp  et  al.,  2005;  Shulman  et  al.,  2005).  The  top  nd  bottom 
panels  of  Fig.  6  are  a  remarkably  good  match  for  Fig.  5  and  L  and 
show  a  more  pronounced  plume  of  cold,  salty  wa  er  moving 
offshore  off  Ano  Nuevo  during  the  August  15  event  than  during 
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Fig.  S.  (A-F)Time  series  of  sea-surface  temperature  (SST.  "C)and  surface  winds  (m/s)  as  sensed  from  the  aircraft  during  August  4  13,  2003.  The  flight  path  s  indicated  by  the 
heavy  black  lines  in  panels  A  and  B  and  the  arrows  indicate  the  magnitude  and  direction  of  the  wind  speed.  (G-L)  As  in  Fig.  5(A-F)  except  for  August  15  "ptember  4,  2003. 


the  August  29  event.  The  model  fields  show  the  Monterey  Bay 
eddy  to  be  much  closer  to  shore  during  the  latter  event,  centered 
near  36.6  N,  122.5  W  with  poleward  flow  on  the  offshore  side. 
These  eddy  circulations  were  responsible  for  the  offshore  advection 
of  the  cold,  salty  plumes  off  Point  Pinos  and  Point  Sur,  and  for 
suppressing  the  offshore  flow  off  Point  Ano  Nuevo  on  August  29. 

While  the  previous  discussion  focuses  on  the  evolution  of  sea- 
surface  temperature,  some  insight  into  the  vertical  structure  of 
the  evolving  upwelling  plume  was  also  obtained  from  the  glider 
formation  flying  experiments  (GFFEs)  described  earlier.  The  first 


GFFE  on  August  6-7  sampled  the  front  when  it  was  young 
and  close  to  the  source  at  Point  Ano  Nuevo  ai  d  the  second 
GFFE  on  August  16  sampled  the  front  when  it  vas  older  and 
had  advected  further  to  the  south  across  the  nouth  of  the 
Monterey  Bay.  In  response  to  the  strong  upwc l ling-favorable 
winds  measured  on  August  6  (Fig.  5C),  three  glider:  were  directed 
that  same  day  to  travel  northwest,  in  an  equilater^  triangle  with 
3  km  length  sides,  towards  the  Point  Ano  Nuevo  re  ,ion  to  sample 
along  the  southwest  side  of  the  upwelling  front  (Fi  ..  7).  Note  that 
this  formation  flew  against  the  prevailing  current  as  sampled  by 
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Fig.  6.  Surface  temperature  (top  panels),  salinity  (lower  panels),and  current  vectors  from  the  NCOM/1CON  model  for  August  15  (left  panels)  and  August  29. 200'  (right  panels). 


all  the  gliders  (blue  arrows  in  the  figure),  which  would  not  be 
possible  using  conventional  profiling  floats  or  drifters.  The  tem¬ 
perature  gradient  estimate  along  the  path  provides  a  higher 
resolution  complement  to  the  SST  (Fig.  5 C)  and  in  the  corresponding 
AVHRR  plots  (not  shown).  The  temperature  gradients  had  the  same 
sign  at  both  depths  shown  (5  and  30  m),  with  warmer  water 
offshore  and  approximately  1  X  cooler  water  inshore.  The  gradient 
was  slightly  stronger  at  30  m.  below  the  surface-mixed  layer  where 
spatial  homogeneity  tends  to  weaken  the  horizontal  gradients. 

During  the  second  GFFE  on  August  16-17  the  gliders  traveled 
from  near  the  center  of  the  bay  towards  the  southeast  for  about 
12  h  and  then  turned  east  and  then  northeast  fora  few  more  hours 
(Fig.  8).  In  this  case  the  gliders  moved  mostly  with  the  prevailing 
current  (blue  arrows)  but  also  across  and  slightly  against  the  flow 
towards  the  end.  The  temperature  gradients  in  the  horizontal 
plane  were  calculated  from  the  glider-observed  temperatures  at  5 
and  30  m  depth  along  the  path  as  in  the  August  6  experiment. 
During  the  first  12  h  the  glider  triangle  was  straddling  the  inshore 
boundary  of  the  upwelling  front,  i.e.,  the  band  of  minimum  SST 
near  122  06'W  (Fig.  5G).  At  5  m  depth,  the  gliders)  furthest 
offshore  measured  cool  temperatures  of  11.5-12  C  in  the 
extended  plume  while  the  gliders  to  the  northwest  measured 
warmer  temperatures  around  13  C  (Fig.  8).  At  30m,  the 
temperature  was  nearly  uniform  (10  C,  dark  blue)  over  much 
of  the  glider  formation  path,  but  with  a  weak  gradient  towards 
the  northwest,  opposite  the  surface  layer.  This  indicates  that  the 
upwelling  filament  was  quite  shallow,  less  than  30  m  deep,  after  it 
had  advected  southward  across  the  Monterey  Bay  By  45  m  depth 
(not  shown)  the  temperature  gradient  was  very  small. 

4  3.  Offshore  structure 

A  central  AOSN-11  issue  is  how  well  the  observing  array  can 
describe  mesoscale  ocean  structure  over  the  region  of  interest. 


Data  assimilating  models  can  be  powerful  tools  for  s  /nthesizing 
the  scattered  and  diverse  AOSN-11  observations  into  a  coherent 
picture,  but  these  models  can  only  describe  that  range  of 
phenomena  that  approximately  obey  the  model  equa  ions  at  the 
resolution  they  are  solved.  Little  is  known  about  wf  it  observa 
tions  are  needed  to  support  an  accurate  and  coherei  t  synthesis 
of  mesoscale  coastal  structure  or  how  these  d  ita  should 
be  statistically  characterized  for  the  assimilation  p  ocess.  The 
observing  power  of  the  AOSN-11  array  is  describi  d  by  data 
collected  over  two  separate  one-week  periods  under  rontrasting 
wind  forcing  (Figs.  9-11).  In  the  first  week  the  offsho  e  sampling 
density  was  not  yet  very  uniform  since  the  glider  were  just 
getting  on  their  assigned  tracks  (Fig.  9A).  During  the  s  cond  week 
some  gliders  were  moved  off  their  assigned  tracks  iy  currents 
(Fig.  QB). 

Winds  over  the  last  week  of  July  and  through  August  6  were 
weak  and  variable  (Fig.  2).  While  the  moorings  pi  avide  time 
series,  the  other  platforms  provide  instantaneous  sai  iples  along 
tracks  and  it  takes  about  a  week  for  them  to  map  out  the  region. 
Even  without  upwelling-favorable  winds  for  a  week  or  nore,  near- 
surface  temperatures  are  2  C  or  more  colder  near  th(  coast  than 
50-100  km  offshore.  The  mean  alongshore-  and  tin  e-averaged 
across-shore  temperature  section  for  this  time  perio  l  (Fig.  10A) 
shows  the  near-surface  temperature  gradient  in  Fit .  9A  to  be 
associated  with  isotherm  slopes  (and  isopycnal  slopes  iot  shown) 
in  the  upper  50-100 m  that  support  geostrophic  ;hear  with 
equatorward  currents  increasing  near  the  surface.  Bel  >w  this,  the 
isotherm  and  isopycnal  slopes  are  reversed,  producing  geos¬ 
trophic  shear  with  poleward  flow  increasing  above  1(  0  m. 

Upwelling-favorable  winds  began  August  6  an  persisted 
through  August  19  (Fig.  2).  Compared  with  the  pre(  -ding  calm 
period,  near  surface  temperature  cooled  about  2  Cv\  thin  10  km 
of  the  coast  as  the  surface  layer  was  blown  offshore  (F  g.  9B).  This 
is  also  evident  in  the  up-warping  of  the  isotherms  within  40  km  of 
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Fig.  7.  Three  gliders  moving  northwest  in  triangular  formation  {gray  lines)  from  18:00:002  August  6  into  late  morning  August  7.  2003.  Colored  lots  indicate  the 
temperature  in  C  at  the  5m  (top)  and  30m  (bottom)  levels;  black  circles  and  lines  show  initial  position  and  formation;  and  open  circles  with  a  cro  s  inside  triangles 
illustrate  the  path  of  the  center  of  the  triangle.  Red  arrows  correspond  to  estimates  of  down-gradient  temperature  in  the  horizontal  plane  along  the  path  c  the  center  of  the 
triangular  formation  The  color  contours  indicate  SST  as  sampled  by  the  TWIN  OTTER  aircraft.  The  blue  arrows  correspond  to  currents  estimated  with  >bjective  analysis 
using  the  gliders’  depth-averaged  current  measurements. 


the  coast  (Fig.  10B).  In  contrast,  the  temperature  section  shows 
that  at  distances  of  50  km  or  more  from  shore  the  near-surface 
temperature  increased  up  to  1  C  as  a  warm  surface  layer  seems 
to  have  accumulated  in  a  weakly  stratified  deep  mixed  layer 
(Fig.  10B).  Below  100m,  no  large-scale  change  in  temperature  or 
temperature  gradient  can  be  detected.  The  accumulation  of  warm 
water  offshore  under  upwelling-favorable  winds  is  surprising 
since  there  is  no  evidence  that  the  wind  strength  or,  consequently, 
the  Ekman  transport  decreases  offshore  or  that  surface  warming 
from  above  increases. 


The  offshore  Spray  gliders  sample  depth-averag  *d  currents  to 
400  m  and  show  a  general,  if  structured,  poleward  fl  >w  offshore  in 
both  time  periods  (Fig.  9).  This  manifestation  of  the  California 
Undercurrent  was  strongest  during  the  weak-wind  eriod  when  it 
took  the  form  of  a  20-km  wide  jet  that  meander*  j  from  25  km 
offshore  crossing  the  southernmost  glider  line  to  (  5  km  offshore 
on  the  northern  line  where  it  exhibited  the  strong  st  currents  of 
nearly  30 cm/s  (Fig.  9A).  Equatorward  wind  in  the  .econd  period 
significantly  reduced  the  poleward  signatures  in  d  pth-averaged 
velocity  (Fig.  9B).  Using  depth  and  alongshore-averaged  current 
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Fig.  8.  Three  gliders  moving  southeast  in  triangular  formation  (gray  lines)  during  August  16  17.  2003.  Colored  dots  indicate  the  temperature  in  Cat  the  S  m  top)  and  30  m 
(bottom)  levels;  black  circles  and  lines  show  initial  position  and  formation,  and  open  circles  with  a  cross  inside  triangles  illustrate  the  path  of  the  center  of  t  ie  triangle.  Red 
arrows  correspond  to  estimates  of  down-gradient  temperature  in  the  horizontal  plane  along  the  path  of  the  center  of  the  triangular  formation.  The  color  co  tours  indicate 
SST  as  sampled  by  the  TWIN  OTTER  aircraft.  The  blue  arrows  correspond  to  currents  estimated  with  objective  analysis  using  the  gliders*  depth  av  raged  current 
measurements. 


as  a  function  of  offshore  distance  to  reference  the  alongshore 
geostrophic  shear  computed  from  alongshore-  and  time-averaged 
sections  of  temperature  (Fig.  10)  and  salinity  allows  cross-sections 
of  absolute  alongshore  velocity  to  be  produced.  These  sections 
show  that  the  decrease  in  poleward  depth-average  current  from 
the  first  plotted  week  to  the  second  resulted  from  offshore 
changes  in  geostrophic  shear  of  alongshore  current  in  the  upper 
150  m  (Fig.  11).  Also  seen  in  Fig.  11  is  how  the  equatorward 
upwelling  jet  in  the  second  time  period  is  largely  confined  to  the 
20  km  closest  to  shore.  More  surprisingly,  aside  from  these 
features,  the  general  pattern  of  alongshore  flow  is  moderately 
complex  but  remarkably  similar  between  the  two  time  periods 
even  though  no  observation  contributes  to  both  time  averages.  It 


remains  to  be  seen  if  this  structure  persisted  for  m  >re  than  16 
days. 

Beyond  the  alongshore-  and  time-averaged  patterns  described 
in  Figs.  10  and  11,  Fig.  9  shows  complex  and  energf  ie  smaller- 
scale  features  in  both  near-surface  temperature  md  depth 
integrated  currents.  Some  of  the  complexity  derives  fiom  plotting 
a  full  week  of  data  together  while  other  features  may  be  a 
consequence  of  internal  tides  that  are  substantially  aliased  by 
shipboard  and  glider  sampling.  But  some  of  the  s  nailer-scale 
features  have  the  appearance  of  eddies  and  fror  ts.  Notable 
are  oppositely  rotating  eddies  found  near  the  mouth  f  Monterey 
Bay.  During  the  first  week  (Fig.  9A)  the  eddy  is  cyclonic  with  a 
clear  poleward  flow  along  the  coast  of  Monterey  Bay  while  it  is 
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Fig.  9.  Raw  observations  of  5-m  temperature  (colored  squares,  scale  to  the  right)  from  AOSN-II  ship  surveys,  moorings  and  gliders  and  depth-averaged  velocity  from  gliders 
(black  vectors,  scale  at  upper  right)  from  (A)  July  31  to  August  6,  2006  and  (B)  August  8-15,  2006.  Spray  gliders  on  lines  extending  well  offshore  sampl  to  the  bottom  or 
400  m  while  Slocum  gliders  on  nearshore  generally  closed  tracks  sample  to  the  bottom  or  200  m. 


anti-cyclonic  in  the  August  8-15  upwelling  period  (Fig.  9B).  The 
anti-cyclonic  eddy  is  implicated  as  a  cause  of  the  spread  of  cold 
water  away  from  Pt.  Pinos  (the  southern  end  of  the  Bay  mouth) 
seen  in  aircraft  SST  as  shown  in  Fig.  5F-H  spanning  August  13-20. 

The  ability  of  the  AOSN-ll  sampling  array  to  describe  the 
smaller-scale  features  hinted  at  in  Fig.  9  depends  on  the  spatial 
and  temporal  scales  of  those  features.  Roughly  put,  objective 


mapping  (Bretherton  et  al.,  1975)  shows  that  mapp  ig  accuracy  in 
the  absence  of  assimilation  into  a  dynamical  moc  *1  depends  on 
the  number  of  observations  within  sampling  volumes  whose 
dimension  is  set  by  the  field's  space  and  time  seal  's.  If  there  are 
enough  points  in  each  such  volume  to  suppress  sit  all-scale  noise 
the  field  can  be  mapped.  In  planning  the  AOSN-1)  field  trial  the 
scales  of  coastal  variability  were  not  well  known.  Fi.  ..  12  shows  the 
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Fig.  10,  Alongshore  and  time-averaged  temperature  measurements  (  C)  averaged 
in  10-m  depth  bins  and  4-km  bins  of  offshore  distance  (shortest  distance  to  the 
coast)  for  (A)  July  31  though  August  6.  2006  and  (B)  August  8-15,  2006.  The  view  is 
southward  with  the  coast  on  the  left 

temporally  lagged  correlation  of  nearly  co-located  temperature 
measurements  and  the  horizontally  isotropic  spatially  lagged 
correlation  of  nearly  simultaneous  measurements  at  the  same 
depth  based  on  all  glider,  shipboard,  and  mooring  data  from 
AOSN-U. 

Several  aspects  of  the  measured  scale  structure  deserve 
comment.  First,  the  spatial  and  temporal  1 /e-folding  scales  of 
15  km  and  2  days,  respectively,  are  nearly  a  factor  of  10  smaller 
than  the  scales  of  typical  deep-water  mesoscale  eddies,  which  are 
commonly  100  km  and  30  days.  Second,  the  correlation  values 
appear  to  increase  rapidly  as  both  the  spatial  and  temporal 
separations  become  small.  Indeed,  at  the  smallest  separations 
of  0.5  days  and  5  km  the  correlation  is  about  0.8,  suggesting  that 
20%  of  the  observed  temperature  variance  has  smaller  scales  than 
this.  Third,  the  correlation  shapes  are  approximately  the  same  for 
all  depths  shown  (and  all  others  that  are  not  shown).  A  more 
complete  discussion  of  the  statistics  of  the  AOSN  11  variability  is 
given  by  Davis  et  al.  (2009)  but  the  main  point  is  that  the  results 
in  Fig.  12  show  that  the  sampling  array  for  the  five  offshore  Spray 
gliders  was  too  thin  to  map  any  but  the  largest-scale  components 
of  the  observed  fields. 

4.4.  Multi-scale  regional  dynamics 

The  multi-scale  dynamical  interaction  between  the  large-scale 
and  mesoscale  variabilities  of  the  August  2003  circulation  and 
their  relationship  with  the  upwelling-related  events  have  been 
investigated  using  the  methodology  multi-scale  energy  and 
vorticity  analysis  (MS-EVA),  and  the  MS-EVA-based  localized 


Fig.  11.  Absolute  alongshore  velocity  (cm/s)  averaged  over  (A)  Ji  y  31  through 
August  6.  2006  and  (B)  August  8  15,  2006.  Geostrophic  shear  Tom  T  and  S 
sections  like  Fig.  10)  was  referenced  by  the  alongshore-averaged  g  der-measured 
depth-averaged  velocity  (from  Fig.  9).  Since  this  constrains  the  ven  cal  integral  of 
velocity.  400m  is  not  a  reference  level  and  the  velocities  there  af  not  weak  by 
construction.  The  velocity  scale  is  at  right  with  positive  numbers  poleward  and 
negative  equatorward.  The  view  is  southward  with  the  coast  on  th  left. 


instability  theory  developed  by  Liang  and  Robinson  (I  305,  2007). 
MS-EVA  is  built  on  the  basis  of  the  multi-scale  windov  transform, 
a  recently  established  mathematical  apparatus  by  Liang  and 
Anderson  (2007).  In  this  framework,  the  multi-scale  pi  nresses  in  a 
given  functional  space  are  represented  on  its  sea. ■?  windows, 
or  mutually  orthogonal  sub-spaces,  each  with  an  exclusive  range 
of  time  scales.  The  objective  of  the  MS-EVA  is,  acc  *rdingly,  to 
investigate  how  energy  and  enstrophy  are  transpo  ted  within 
their  individual  scale  windows,  redistributed  throug  the  inter 
action  between  the  windows,  and  converted  and  dissi  >ated  when 
a  system  is  steered  forward.  The  redistribution  is  characterized  by 
a  concept  called  perfect  transfer ,  which  has  been  rigorously 
connected  to  the  two  important  processes,  barotropi-  instability 
and  baroclinic  instability,  in  geophysical  fluid  dynami  s. 

The  AOSN-I1  Harvard  Ocean  Prediction  Syst(  m  (HOPS) 
model  fields  have  been  analyzed  in  three  mutually  ext  usive  time 
sub-spaces:  a  large-scale  window  ( >  8  days),  a  mesosc  tie  window 
(0.5-8  days),  and  a  sub-mesoscale  window  (<0.5  days).  The 
analysis  has  been  performed  during  the  upwelling-  elaxation- 
upwelling  cycle  of  winds  shown  in  Fig.  2.  Upwellir  ^-favorable 
winds  occur  from  6-19  August  to  23-31  August  and  th  relaxation 
period  is  19-22  August.  The  dynamics  of  the  fields  are  driven 
externally  by  wind  and  surface  fluxes  and  are  cl  aracterized 
internally  by  the  transfers  of  available  potential  energ  /  (APE)  and 
kinetic  energy  (KE)  from  the  large  scale  to  the  mesost  ale  both  (i) 
during  upwelling  and  (ii)  during  relaxation.  Off  of  P(  nt  Sur  the 
wind  destabilizes  the  system  directly  during  upwelling .  with  wind 
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Fig.  12.  Temporally  lagged  correlation  (left)  and  spatially  lagged  isotropic  correlation  (right)  of  temperature  at  various  depths.  Temporally  lagged  result  ir  eludes  all  AOSN-U 
observations  that  are  within  5  km  of  each  other  while  spatially  lagged  correlation  includes  all  data  whose  observation  time  differs  by  less  than  12  h.  The  si  nallest  lag  bins  do 
not  include  data  correlated  with  themselves.  The  legend  at  upper  right  describes  the  color  and  symbol  code  for  different  depths  in  meters. 
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Fig.  13.  Energy  transfer  term  from  the  HOPS  large-scale  window  to  mesoscale  window  at  10  m  Top  row  available  potential  energy  (APE);  bottom  row  k  netic  energy  (KE). 


energy  captured  within  the  large-scale  window  and  continually 
released  to  the  mesoscale.  In  contrast,  near  and  in  Monterey  Bay. 
the  wind  tends  to  stabilize  the  southward  coastal  current,  energy 
is  stored  in  the  large  scale,  and  instability  occurs  only  when  the 
external  wind  constraint  is  relaxed. 

Fig.  13  illustrates  the  available  potential  and  kinetic  energy 
transfer  from  the  large-scale  window  to  the  mesoscale  window  at 
4  day  intervals  from  13  to  25  August  and  Fig.  14  schematizes  the 


energy  transfer  process.  During  (i)  upwelling  (1  .17.  and  25 
August)  the  greatest  direct  transfer  of  APE  and  KE  i<  in  the  area  to 
the  west  of  Point  Sur.  Once  (ii)  relaxation  occur.  (21  August), 
energy  transfer  in  this  location  is  significantly  ret  uced  and  the 
location  of  energy  transfer  moves  to  the  region  north  of  Monterey 
Bay.  near  Point  Ano  Nuevo.  The  generated  mesost  lie  structures 
propagate  northward  along  the  coast  with  dispers  on  properties 
similar  to  those  of  free  thermocline-trapped,  coastally  trapped 
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(A) 
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Fig.  14.  Schematic  of  the  major  dynamical  processes  in  the  upper  layer  near  Monterey  Bay  during  (A)  upwelling-favorable  winds  and  (B)  relaxation.  Major  energy  transfer 
from  the  large-scale  window  to  the  mesoscale  window  takes  place  west  of  Point  Sur  during  upwelling.  while,  during  relaxation,  major  energy  transfer  take  place  north  of 
and  in  Monterey  Bay.  The  generated  mesoscale  events  propagate  northward  in  the  form  of  coastally  trapped  waves. 
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Fig.  15.  Four  snapshots  of  particle  dispersion  in  the  Monterey  Bay  in  the  vicinity  of  an  LCS  (red  line)  during  upwelling  favorable  winds  The  LCS  and  particle  <  ispersion  were 
computed  from  pseudo-drifters  inserted  into  the  CODAR  HF  radar  surface  current  maps  obtained  during  the  AOSN-11  experiment  Results  are  shown  for  (A)  /  jgust  13.  2003 
at  0700;  (B)  August  14,  2003  at  0100:  (C)  August  14.  2003  at  2100,  and  (D)  August  15.  2003  at  0700.  For  additional  discussion  of  this  figure,  see  (Shadde  i  et  al..  2009). 


waves.  These  waves  have  been  previously  observed  off  Point  Sur 
by  other  investigators  (Ramp  et  al.,  1997). 

4.5.  Lagrartgian  coherent  structures  in  the  Monterey  Boy 

Lagrangian  Coherent  Structures,  or  LCS,  is  a  computational 
technique  for  processing  velocity  field  data  in  a  flow,  including 


currents  in  the  ocean  (Shadden  et  al.,  2005).  The  basi  methodol 
ogy  for  this  technique,  based  on  the  idea  of  finite  tir  ie  Liapunov 
exponents  (measuring  how  particles  are  separating  in  forward 
time,  or  converging  in  backward  time)  originated  du  ing  the  few 
years  before  AOSN-U  (Haller.  2002;  Lekien,  200-  )  and  was 
subsequently  further  developed  using  the  AOSN-11  da>a. 

One  of  the  things  LCS  can  provide  is  to  disceri  regions  of 

different  particle  (or  drifter)  behavior.  In  conjunction  with  model 
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predictions,  this  can  be  a  useful  predictive  tool.  The  interesting 
fact  is  that  even  in  a  complex  flow  such  as  that  in  Monterey  Bay, 
LCS  provides  time-varying  curves  that  separate  particles  that  have 
different  dynamical  fates  and  particles  do  not  cross  these  LCS. 
There  is,  for  instance,  a  well-defined  curve  (technically,  a  repelling 
LCS)  that  separates  particles  in  the  inner  bay  from  those  in  the 
outer  bay  (Fig.  15).  Even  though  this  curve  is  moving  in  time,  and 
even  though  the  flow  is  complex  with  multiple  space  and  time 
scales,  this  curve  on  the  scale  of  the  bay  itself  sharply  separates 
particles  that  have  different  fates.  Particles  on  one  side  of  the 
curve  stay  on  that  side  Thus,  the  fluid  to  the  right  of  the  red  curve 
may  be  regarded,  for  this  time  period,  to  be  re-circulating  in  the 
bay.  That  there  even  is  such  a  well-defined  curve  is  remarkable. 
Numerical  evidence  supports  the  idea  that  this  curve  (although 
perhaps  not  LCS  for  smaller-scale  structures)  is  robust  to  data  and 
model  uncertainty. 

Since  the  problem  is  time-varying,  the  recirculation  zone  does 
not  persist  indefinitely,  and  how  such  structures  appear  and 
disappear  is  related  to  the  mechanisms  of  mixing  and  transport  in 
the  ocean.  LCS  can  determine  the  boundaries  of  vortex  regions  and 
elucidate  the  mechanisms  of  entrainment  and  detrainment  from 
vortex  structures.  It  does  so  with  greater  clarity  and  precision 
than,  for  example,  vorticity  plots,  as  has  been  demonstrated  in 


computations  with  laboratory  vortex  rings  (Shadden  et  al„  2006). 
An  example  for  the  ocean  has  been  computed  using  output  from 
the  JPL/ROMS  model  (Fig.  16).  The  fluid  in  the  lobes,  defined  as  the 
regions  initially  bounded  by  the  LCS,  was  initially  distributed  as 
in  Fig.  16B.  Over  time,  the  fluid  moved  away  from  the  repelling 
LCS  (orange)  and  towards  the  attracting  LCS  (blue)  ultimately 
being  entrained  by  the  blue  vorticies  as  shown  in  figs.  16C  and  D. 
Transport  and  mixing  occurs  via  the  motion  of  the  lobes  since  no 
fluid  crosses  the  LCS.  That  drifters  in  the  green  lobes  are  flushed 
out  to  sea  while  drifters  in  the  brown  lobes  recirculate  near 
the  coast  evinces  the  property  that  LCS  separatt  regions  with 
different  particle  fates  and  residence  times.  This  mnkes  LCS  a  very 
promising  tool  for  studying  the  distribution  of  water-mass 
properties,  determining  the  fate  of  pollutants  and  oil  spills,  and 
planning  how  to  seed  surface  drifters  and  autonomous  vehicles 
for  field  studies  of  ocean  dynamics. 

4.6.  Ocean -leaving  radiance 

The  ocean-leaving  radiance  data  were  sometii  les  helpful  in 
delineating  the  water-mass  boundaries  in  addition  o  showing  the 
distribution  of  marine  organisms,  particularly  durii  g  weak  winds 
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Fig.  16.  LCS  along  the  central  California  coast  near  the  Monterey  Bay  as  computed  from  the  JPL/ROMS  model  output  The  LCS  are  depicted  for  {A  and  E  July  25,  2006  at 
0900.  (C)  July  26.  2006  at  0100.  and  (D)  July  26.  2006  at  2300.  Intersections  of  the  attracting  LCS  (blue  curves)  and  the  repelling  LCS  (orange  curves)  define  he  boundaries  of 
lobes  colored  green  and  brown. 
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Fig.  17.  Chlorophyll-a  distributions  as  sensed  from  the  hyperspectral  sensor  on  board  the  TWIN  OTTF.R  aircraft  for  (A)  August  IS.  (B)  August  20.  and  (C)  A  igusi  29,  2003. 


when  a  thin  layer  of  warm  surface  water  masked  the  underlying 
structures.  The  hyperspectral  data  were  collected  using  a  HOBI 
Labs  HydroRad  hyperspectral  radiometer  measuring  downwelling 
irradiance  (Ed),  upwelling  radiance  (lu),  and  downwelling  radi¬ 
ance  (Ld).  From  these  parameters,  the  water-leaving  radiance  (Lw) 
was  calculated  using  different  algorithms  for  a  cloud-covered  or 
cloud  free  sky  as  appropriate.  The  water-leaving  radiance  from 
four  spectral  channels  was  combined  subsequently  to  estimate 
the  chlorophyll-a  (Chl-a)  content  of  the  surface  water  using  the 
SeaWiFS  OC4  version  4  algorithm.  Details  of  the  calculations  can 
be  found  in  a  separate  technical  report  (Anderson  and  Ramp, 
2005). 

Three  examples  of  the  surface  distributions  of  Chl-a  from  15, 
20,  and  29  August  were  sampled  during  the  first  upwelling  event, 
the  wind  relaxation,  and  the  second  upwelling  event  (Fig.  17). 
Comparison  with  Figs.  5G,  H,  and  L  shows  almost  no  Chl-a 
seaward  of  the  17  C  isotherm,  which  may  be  taken  as  the 
boundary  between  the  California  Current  water  and  the  coastal 
water.  The  exception  was  an  offshore  patch  of  high  chlorophyll 
near  36.9  N,  123.0  W  on  29  August,  which  according  to  SST 
(see  Fig.  5K)  was  located  in  old  upwelled  water  which  had  been 
previously  advected  offshore  on  25  August.  All  the  Chl-a  maps 
were  patchier  than  SST  with  the  highest  values  occurring  away 
from  the  coast  in  older  upwelled  water.  The  vigorous  advection 
and  mixing  during  active  upwelling  events  apparently  prohibited 
the  static  stability  required  for  a  phytoplankton  bloom  to  take 
place.  Once  the  wind  died  down  and  surface  warming  re-stratified 
the  water  column,  growth  accelerated  rapidly.  The  offshore 
“squirt”  south  of  the  bay  (Fig.  5H)  contained  patches  of  high 
Chl-a  water  being  advected  offshore  from  Point  Pinos  by  the  eddy 
identified  by  the  gliders.  The  Chl-a  distributions  on  15  and  29 
August  show  the  greater  offshore  advection  of  material  during  the 
earlier  event  vs.  the  later  one.  The  nearly  oligotrophic  water  close 
to  shore  off  Ano  Nuevo  on  29  August  (dark  blue  in  Fig.  17C)  further 
supports  the  notion  that  the  California  Current  confined  the 
upwelling  center  to  the  nearshore  region  during  the  second  event 
and  encouraged  the  southward,  rather  than  offshore  advection 
of  the  cold,  nutrient-rich  water. 


5.  Impacts  and  conclusions 

The  AOSN-II  experiment  successfully  demonstrated  the  viabi¬ 
lity  of  several  new  technologies  and  techniques  in  addition  to 
illuminating  some  basic  scientific  issues  surrounding  coastal 
upwelling  off  central  California.  Gliders  running  closed  cells  and 
racetracks  were  able  to  determine  robust  statistics  of  the  flow  due 


to  their  impressive  persistence  and  spatial  coverag  .  They  also 
observed  smaller-scale  features  such  as  the  eddy  at  t  le  mouth  of 
the  Monterey  Bay.  Gliders  flying  formations  were  ab  e  to  sample 
frontal  gradients  in  both  the  horizontal  and  vertical  w  hile  moving 
across  or  against  the  flow  field  as  well  as  moving  w  th  it  These 
observations  were  not  previously  possible  with  float  or  drifters, 
and  were  not  possible  without  an  impressive  level  of  utonomous 
control  over  the  small  glider  fleets.  The  low-flying  lircraft  was 
able  to  map  out  the  three-dimensional  structure  of  th  wind  fields 
with  a  level  of  detail  sufficient  to  make  quantitative  ■  omparisons 
with  the  latest  high-resolution  atmospheric  mode! ,.  By  flying 
beneath  the  clouds,  the  aircraft  also  collected  a  tine  series  of 
synoptic  snapshots  of  the  surface  temperature  and  o<  ean  leaving 
radiance  fields  that  were  not  available  any  other  way.  The  SST  was 
a  key  parameter  for  assimilation  into  all  three  nume  ical  models 
on  a  daily  basis.  The  Lagrangian  Coherent  Structures  (I  CS)  analysis 
technique,  relatively  new  to  the  field  of  oceanography,  showed 
great  promise  for  understanding  the  flow  dynamics  md  particle 
fates,  whether  using  observational  data  (surface  ci  rrent  maps 
from  HF  radar)  or  numerical  model  (JPL/ROMS)  outp  t. 

Using  the  combined  power  of  these  new  methoc  ologies,  the 
following  picture  of  the  circulation  in  and  around  tl  e  Monterey 
Bay  during  a  series  of  upwelling  and  relaxation  events  has 
emerged.  Two  upwelling  events  were  observed,  on  lasting  14 
days  and  the  second  nine,  with  a  3-day  relaxatu  n  event  in 
between.  The  wind  speed  was  similar  between  the  t\A  d  upwelling 
events,  blowing  steadily  from  the  northwest  at  10-15  n/s.  During 
the  relaxation  winds  were  poleward  at  about  5  m;.  All  three 
events  were  well  sampled  by  the  aircraft,  gliders,  and  ships. 

Both  upwelling  events  began  when  a  pool  of  coh  water  less 
than  13  C  appeared  at  the  surface  off  Cape  Ano  Ni  evo.  As  the 
wind  continued  to  blow  for  about  5  days,  the  cold  po  >1  expanded 
and  spread  offshore  and  to  the  south  across  the  nr outh  of  the 
Monterey  Bay.  A  key  difference  between  the  two  events  was 
the  degree  of  offshore  spreading;  The  first  event  spr  ad  offshore 
and  to  the  south  while  the  second  event  spread  exclu  ively  to  the 
south.  Both  models  and  observations  indicate  tl*  it  the  key 
difference  between  these  events  was  caused  by  the  position  of 
the  offshore  meanders  and  eddies  of  the  Califor  lia  Current 
System.  During  the  intervening  relaxation  event,  tf  s  boundary 
between  the  coastal  and  CCS  waters,  as  indicated  I  y  the  17  C 
isotherm  and  the  Chl-a  content,  moved  40  km  onsho  e  to  within 
5  km  of  the  coast  at  its  closest  incursion.  This  f  ontal  zone 
effectively  blocked  the  offshore  advection  of  cold  watt  during  the 
second  upwelling  event,  and  the  equatorward  geosi  ophic  flow 
along  the  front  aided  the  southward  advection  of  th<  cold  water 
by  the  wind.  The  plume  thinned  vertically  as  it  moved  southward 


S.R.  Ramp  er  a/.  /  Deep-Sea  Research  II  56  (2009)  68-86 


B5 


across  the  bay  to  less  than  30  m  at  the  southern  side.  The 
dynamics  of  the  eddy/frontal  movements  remains  a  mystery  as 
they  were  not  always  correlated  with  the  local  wind  stress.  While 
the  front  clearly  moved  onshore  during  the  relaxation,  it  also 
moved  onshore  at  other  times  and  did  not  retreat  offshore  during 
the  second  upwelling  event.  Using  the  entire  glider  data  set,  the 
temporally  and  spatially  lagged  data  indicate  e-folding  scales  of  2 
days  and  15  km,  respectively  for  these  eddies,  much  smaller  than 
the  scales  (30  days  and  100  km)  typically  observed  in  deep-water 
eddies  offshore.  The  California  Undercurrent  was  present  in 
approximately  the  same  location  throughout  the  experiment, 
further  offshore  in  the  north  than  the  south.  The  changes  in 
the  vertically  averaged  currents  were  due  primarily  to  changes  in 
the  upper  water  column  rather  than  changes  in  the  CUC. 

These  results  suggest  directions  for  future  dynamical  studies. 
The  spatial  distribution  of  the  wind-stress  curl  suggests  that 
positive  curl  on  the  inshore  side  of  a  sub-mesocale  atmospheric 
jet  may  be  important  in  driving  upwelling  in  the  Monterey  Bay,  in 
addition  to  divergence  at  the  coast.  The  HOPS  model  suggests  that 
the  energy  transfer  from  the  large  scale  to  the  mesoscale  takes 
place  via  mixed  barotropic/baroclinic  instability  off  Point  Sur 
during  upwelling  events  and  in  the  northern  Monterey  Bay  during 
relaxations.  More  focused  observational  studies  are  needed  to 
verify  these  ideas. 
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